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Abstract 

The assisted tunneling of a metastable state between barriers is investigated analyt- 
ically by means of a simplified one dimensional model. A time dependent perturbation 
changes the pole spectrum of the wave function introducing a larger decay constant. 
New insights about the decay of a metastable state are found. The scheme is exempli- 
fied for parameters corresponding to the nuclear process of a decay. 
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1 Introduction 



Quantum tunneling is the conventional theoretical paradigm for the explanation of a gamut 
of atomic, molecular and nuclear phenomena. Classically forbidden regions can be accessed 
by a quantum object whose behavior is described in terms of a wave function. 

Shortly after the advent of quantum mechanics, Gurney and Condon and simultaneously 
Gamow, explained the huge differences between a decay lifetimes of similar unstable nuclei 
using the concept of tunneling through a barrier. [1, 2] This framework has remained as 
reliable today as when it was proposed almost a century ago. [3, 4] The question of assisting 
tunneling, i.e. accelerating or decelerating the decay process is of the utmost importance 
in the nuclear context and elsewhere. In the nuclear case, acceleration of the decay process 
could change dramatically the treatment of radioactive waste as well as providing alternative 
fuel sources. 

Analytical expressions for the various tunneling processes are extremely important. It is 
virtually impossible to follow numerically the evolution of a wave function to the long times 
involved in decay processes. There is a huge gap between the natural time scale of nuclear 
phenomena of the order of and the decay time scale of the order of milliseconds 

to millions of year. [6] The situation is somewhat, but not radically different, for atomic 
phenomena. 

Not surprisingly, assisted tunneling is extensively investigated for solid state systems. 
The main theoretical tool in this endeavor is the Floquet formalism[7]. This method in- 
volves large matrices and it is computationally quite intensive. [8]. Other methods such as 
the "elevator effect'' model of resonance assisted activation attempt to attack the problem 
in a semi-analytic manner [10]. Ivlev[9] has recently developed a complex time method to 
investigate assisted tunneling, aimed at a nonperturbative approximate treatment of the 
tunneling through special barriers. 

There seems to be a clear consensus, that exposing a system to external excitations, time 
harmonic or not, can enhance its tunneling rate. The absence of exact analytical results, 
and reliable long term numerical results, especially for the nuclear case, makes difficult the 
assessment of experimental feasibility of implementing assisted activation. 

In the present work we investigate an extremely simple model of a one dimensional quan- 
tum mechanical system located initially between fixed barriers. A time harmonic potential 
such as a low frequency -as compared to the natural frequency of the system- electric field, 
is then applied. The simplicity of the model allows a full, albeit perturbative, exact ana- 
lytical solution. The perturbation causes a qualitative change in the pole structure of the 
wave function. The spectrum of poles arising from the normalization of the wave function 
that rule the decaying behavior of the metastable state, is extended. These poles in the 
complex momentum space, appear as the wave function is expanded in a complete set of the 
unperturbed system. The poles contribute in a contour integration of the wave amplitudes 
in momentum space. [4]. They are not introduced ad-hoc, as Gamow did with his complex 
energy method, but, arise in a systematic expansion of the wave function naturally. The 
time development of the wave function is dominated by the lowest energy poles. The intro- 
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duction of the perturbation opens a new line of poles. Here again only one of them appears 
important. This new pole produces a bigger decay constant. 

The result is not specific to time varying fields. The new line of poles arises even for 
a static spatially inhomogeneous potential such as that of a constant electric field. In a 
nuclear setting, such a static field will be certainly almost completely screened by the atomic 
electrons, whereas for a time varying potential, the efficiency of such screening is diminished. 
A discussion of the electronic screening problem is beyond the scope of the present work. 
We henceforth treat only a time harmonic perturbation that vanishes in the zero frequency 
limit. 

The external vibration shakes the system back and forth between barriers and accelerates 
the tunneling process. We analyze the general and, mostly qualitative results, based on the 
pole spectrum. New insights into some troubling aspects of metastable state decay are offered 
in section 5. Extensive mathematical details will be presented in a forthcoming work. 



2 The Model 

The Schrodinger equation for the one dimensional system located in between barriers is taken 
to be 

9* -1 

Comparing to a finite size square barrier, A represents the height multiplied by the width. 
For nuclear a decay A 2, m ^ 20/m~^, xq ^ 10 fm, mX xq = 400 >> 1 

We choose as initial nonstationary wave, a gaussian wave packet. This allows a full 
analytical treatment. In a decay, it would correspond to a preexisting spatially symmetric 
state. 



^i){x,t = Q) = N (2) 

with N a constant factor normalization, and A is the width parameter of the packet. Inside 
the barrier region the packet disperses if A < xq. , but, soon enough it stops dispersing 
due to the presence of the barriers. It can then spread only through the tunneling process 
governed by the barriers. 

The even (e) and odd (o) stationary states of eq.(l) are readily found to be 

(h.\ ( \ — j cos{k x) if \x\ < Xq . s 

ne\ ) Xe[x) - I ^ cos(/i; x) ± Bsin{k x) if x > xq or x < -xq ^ ' 

, ( \ _ { sin{k x) if \x\ < Xq , . 

no[ ') Xo[x) - I _|_^ ^^^^^ _^ Dsin{k x) if x > xq or x < -xq ' ^ 
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The set of even-odd functions is orthonormal and complete. The normahzation factors^ 
are extremely important and determine the location of the poles in the complex plane. As 
shown below, these poles govern the exponential decay of metastable wave functions. 

The normalization factors are 

= Ti {A{kf + B{kf) 
A{k) = 1 - sin{2 k Xq) 



2 k 



B(k) = ^ (cos(k xo)f 
k 



e {ne{k)f = [{k- sin{2 k Xq) + (mA {cos{k Xo)ff^ (5) 

K(fc))2 = 7:{C{kf + D{kf) 



C{k) = {sin{k xo)f 

D{k) = 1 + sm(2 k Xq) ^ 



k^ {no{k)f = TT ((A; + sin{2 k Xo) + (mA {sin{k xq))')') (6) 

The time harmonic potential perturbation reads 

V{x,t) = jj, X sin{uj t) (7) 

with ^ a coupling constant. For an external electric field of intensity £"0 interacting with an 
a particle of charge 2 |ge|, A* = 2 \qg\ Eq. 

The wave function is expanded in the complete set of even and odd states 

^'(a;,t)=V/ Xi{k,x) ai{k,t) dk (8) 

dots denoting derivatives with respect to t. 

The Schrodinger equation for the amplitudes ae,o becomes 

e ^ < Xe{k,x)\ V{x,t)\Xo{k',x) > dk' ao{k',t) 

e" <Xoik,x)\V{x,t)\xeik',x)>dk' ae{k',t) (9) 



^The determination of the normalization and the completeness issue will be dealt with in an extended 
version of the paper. [5] 
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The leading contribution to the matrix element of the interaction can be evaluated exactly 



/oo 
Xe{k,x) V{x,t)xo{k',x) dx 
-oo 



1 d5(k-k') . , , 



Inserting eq.(lO) in eq.(9) we obtain 



i de{k,t) = -— —{a'^{k,t) - n'^{k)/no{k) ao{k,t) ao{k,t))ii sin{u t) 

no[k) ne[k) m 

i do{k, t) = ^ {a'^{k, t) - n'^{k)/ne{k) ae{k, t) - ^-^ae{k, t))^ sin{uj t) (11) 

primes denoting derivatives with respect to k. 
Rescaling to dimensionless variables 

k ^ k^ ^ 



t — > i — CO t 

- ~^=^== (12) 



m 



eq.(ll) becomes 



TT 



i ae{k,t) = -jl — — -~-{a'^{k,t) - n'^{k) / no{k) ao{k,t) - i kt ao{k,t)) sin{t) 

no[k) ne[k) 

i ao{k,i) — —jl = —{a'^{k,t) — n'^{k) / ne{k) ae{k,i) —iki ae{k,i)) sin{i) (13) 

no{k) ne{k) 

primes denoting derivatives with respect to k and dots representing derivatives with respect 
to i. 

jl is the relevant parameter of the problem. It is small for external angular frequen- 
cies w^jjj > ix^ m. ^ For lower frequencies, nonperturbative solutions are required. The 
smaller the electric field, the lower the frequency for which the perturbative solution will be 
appropriate. 

The amplitudes of eq.(13) are further expanded in powers of jl 



^For example, in a decay case with w 1 eV"^ corresponding to a strong electric field amplitude of 



Eo'^b IQ" Volt/m, Umin > 10 sec~ 
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n=0 
oo 

ao{Kt)^Y. (14) 

n=l 

with initial conditions 

l){n){k^t = Qi) = Oforn>l 

bf\k,i) = boik)/neik)^ Ne-^Xe{k)dx (15) 

J—oo 

We extracted the normahzation factor of the even wavefunctions ne{k) for the sake of conve- 
nience. It is straightforward to show that all the odd n powers of ae{k,i) vanish identically, 
and the same is true for the even n powers ao{k, t). The substitution of eqs.(14) into eq.(13) 
produces a separate equation for each order n. The equations are integrable exactly order 
by order, although the expressions of higher order amplitudes become increasingly involved. 



3 Poles of the Wave Function 

Figure 1 shows — ttttt, and — -— , for the parameters mX xq — 400, = 10/m. 

k^ ne{ky k^ np[ky 

The thin and tall spikes are due to the extreme closeness of the minima of the normalization 
factors to their complex zeroes. 

The amplitudes of eq.(14) determine the wave function through eq.(8). Integration over 
k is dominated by the pole structure of the normalization factors rio^eik). 

The normalization factors of eq.(5,6) are even in the argument k. Consequently, the 
behavior around their zeroes can be expressed in the form 



^ 1^J> {k' - k^fr + (16) 

where j enumerates the pole number. 

For the 5 barriers we have chosen, and in the limit of mX xo » 1, we can find simple 
expressions for the locations of the poles as well as the value of the norm factors on the real 



axis at k — k^l- 



(2n + 1) TT mX n TT mX 



2(1 + mX Xq) ° (1 + mX Xq) 

2 mX xKmX xq + 4) mX Xo(mA xq + 4) 

" (271+1)%^ = ^ 

^^^(2n + ir^ = nr) 

16 m^X^ Xq m?X'^ Xq 
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Figure 1: -^-^-j^, and ^2 n (k)'^ ' ^ ^ function of k in units of fm for the parameters 
mX xo — 400, xq — 10 fm. 
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The poles Qn, Qm are located symmetrically above and below the real momentum axis at 
k Xq approximately equal to an odd multiple of — for the even case and, and a multiple of 
TT for the odd case, where the bound states of the infinite wall case lay. 

^ _ / (2n + l)7r mX y . 
^" ~ \2 (l + m\xo)) ' 




^" = ((TT^) 

The imaginary parts of the poles are orders of magnitude smaller than the real parts. 
For example, inserting the parameters corresponding to a decay, the first even pole appears 

at C[e,n=l 

fti 0.024 ±i 3* 10'^ fm^^. For a finite size width, the imaginary parts are much 
smaller. The sharp spikes seen in figure 1 lead to the dominance of the poles in the spectrum 
of momenta. 

When eq.(8) is evaluated by contour integration in the complex k"^ plane, the contour has 
to be closed from below. In the lower half-plane the convergence is insured by the exponential 

e 2 m . 

The poles are separated from each other, their contribution adds up. The negative 
imaginary part of each pole induces a time decaying exponential. Each exponent determines 
a different decay constant and decay time. If the original wave function is even in space, 
as in the expression of eq.(2),only even poles contribute to the unperturbed decay process. 
Due to the initial confinement of the wave to the intcr-barrier region, the most important 
contribution arises from the first even pole. The infiuence of higher order poles is hindered 
by the wave packet transform to momentum space 



^{k) - (19) 

Therefore, the first pole drives the decay process and determines the lifetime of the 
quasistationary state. 

4 Acceleration of the Decay Process 

The structure of the solution of eq.(ll) for the amplitude series of eq.(14) and to the lowest 
order in /i for the even amplitude is found to be 



ae{k, i) = ne(k) bo{k) ( ^ + Z^' i)) 
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The odd amplitude reads 



^ k 64 (^) K + e5{i) < 



+ ^ + ^4^ (20) 



mj) -'^^^ (21) 



where we have extracted a normalization factor in advance of the integration over the even- 
odd set of wave functions x(fc), as we want to identify the influence of the pole structure on 
the amplitudes. In eqs.(20,21) we have suppressed the k dependence of the normalization 
factors in the denominators for brevity. 

We can readily analyze the structure of the perturbed amplitudes. The amplitudes in 
eqs. (21,20) differ from the unperturbed amplitude bo{k) by the appearance of poles for the 
odd set of wave functions introduced by powers of no{k) in the denominators. 

The first even pole of eq.(17) generates a decay constant of the form 

Ae = 

z m V 7i 

(22) 




8 m Xq m? 



whereas the first odd pole generates a decay constant 



Ao = 8 Ae (23) 

The factor of 8 is specific to the 5 model. In general there should be roughly this order of 
magnitude increase in the decay constant. 

At the same time, the real part of the pole that enters the wave function mainly through 
the momentum space transform of the initial wave function of eq.(19) changes also 



e i« < 



(24) 
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The left hand side in both hnes in eq.(24) correspond to the first even pole contribution 
whereas the right hand sides pertain to the first odd pole and the second even pole respec- 
tively. 

There arise then two effects: An enhancement of the decay constant by a factor of 8 
and, a damping factor. A rough estimate of the importance of the news terms can now be 
performed. 

If initially A ^ Xq, then the first odd pole contributes to the wave function about 10%, 
while the second even pole contributes much less than 1%. Even without getting into the 
cumbersome details of the wave function, the survival probability in the inner region between 
barriers will be influenced markedly by the introduction of the perturbation. The wave will 
tunnel faster. 



5 Flux Continuity 

The question of conservation of probability or fiux continuity seems unavoidable at this 
junction. This problem arises even before adding a harmonic perturbation. If the wave 
diminishes with time everywhere it appears as if unitarity is broken. This should not occur 
in the evolution of a wave function with the Schrodinger equation and a real potential. The 
contradiction arises in other conventional approaches to the decay of a metastable state [4]. 
In the scattering method the paradox is apparently resolved by exhibiting another piece in 
the wave function that increases exponentially with time. 

However, in the present treatment of the problem, the wave function receives contribu- 
tions from the poles in the lower momentum half-plane only. There seems to be no time 
increasing piece at all. The exponentials originating from those poles is decreasing in time. 
On the other hand, fiux has to be continuous and overall probability should be conserved. 
The answer to this conundrum brings to the fore the mathematical beauty of the decay 
process. 

Consider the two separate regions, the inside zone around a; ~ and the long distance 

zone of a; — s> cxD, far away from the microscopic location of the barriers. In the inner region 
the integration over momentum space is dominated essentially by the poles. This is true 

until very long times. At very long times the exponential factor e * 2™ oscillates wildly for 

small changes in the momentum. The integral of eq.(8) will start to suppress the importance 

of the pole share in the integral. The pole spike will be twisted by the oscillating factor 

with alternating positive and negative parts. It is possible to estimate the time at which the 

assumption of pole dominance fails as follows. The width of the pole spikes read off from 

eq.(16) is 5{k'^) ~ 7. For /c^ ft; 7 the inverse of the normalization factor essentially vanishes 
1 

^ 0. The pole dominance will end when 



k"^ t -ft 
2 m 2 m 



(25) 
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Inspection of the values of 7 listed in eq.(17), it is found that this time is of the same order 
of the lifetime — . Hence, in the inner region, the wave will decay until a time of the order of 
the lifetime. After that time, other contributions to the wave arising from momenta outside 
the poles have to be included. The wave leaks out of the inner region for very long times of 
the order of the lifetime of the state. 

The behavior in the long distance region is completely different, even for short times. 
For X — > 00 the oscillations in the integrand of eq.(8) stem from the harmonic functions in 
space X of eqs.(3,4). For the long distance region, the poles will never dominate , neither at 
short times nor at asymptotically long times. The wave is built from the whole momentum 
spectrum. Moreover, as time passes and, in the spirit of the stationary phase approximation, 
smaller and smaller momenta will dominate. This effect will enhance the wave at long 
distances, due to the momentum space factors of eq.(19). Consequently, the wave function 
will decay in the inner region and grow in the outer region and probability will be conserved. 
The same is true when more than one decay constant is present. 

6 Conclusions 

The investigation of assisted tunneling carried out in the present work implies that, there 
are distinctive qualitative features induced by the introduction of a time dependent per- 
turbation. The main signal for the process is the appearance of another decay constant, 
bigger than the one of the unperturbed case. The state decays as if it were composed of two 
channels with different partial decay widths. A large component with a longer lifetime and 
a small, but presumably non- negligible component, with a much smaller -eight times smaller 
for mA Xo >> 1- lifetime. Evidently, the fine details are important for an accurate determi- 
nation of the wave function, and will be shown in a coming publication. The reverse process 
of tunneling into a region enclosed between barriers could be affected by the mechanism of 
assisted tunneling also. Such a possibility is of relevance for atomic and nuclear processes, 
such as nuclear fusion. 
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